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Abstract. We study the motion of a particle embedded in a time independent 
periodic potential with broken mirror symmetry and subjected to a Levy noise 
possessing Levy stable probability law (Levy ratchet). We develop analytical ap- 
proach to the problem based on the asymptotic probabilistic method of decompo- 
sition proposed by P. Imkeller and I. Pavlyukevich [J. Phys. A 39, L237 (2006); 
Stoch. Proc. Appl. 116, 611 (2006)]. We derive analytical expressions for the 
quantities characterizing the particle motion, namely the splitting probabilities of 
first escape from a single well, the transition probabilities and the particle current. 
A particular attention is devoted to the interplay between the asymmetry of the 
ratchet potential and the asymmetry (skewness) of the Levy noise. Intensive nu- 
merical simulations demonstrate a good agreement with the analytical predictions 
for sufficiently small intensities of the Levy noise driving the particle. 

1 Introduction 

Levy motion, also referred to as "Levy flights" , stands for a class of non-Gaussian Markovian 
random processes whose stationary increments are distributed according to the Levy stable 
probability laws originally studied by French mathematician Paul Pierre Levy [Tj. The term 
"Levy flights" (LF) was coined by Mandelbrot [2], who thus poeticized this type of random 
motion which is now considered as a paradigm of non-Brownian random walk. Similar to Brow- 
nian motion, LFs have a solid probabilistic background. Indeed, the central limit theorem and 
the properties of Gaussian probability laws and processes constitute mathematical foundation 
of the Brownian motion [3]. At the same time, the generalized central limit theorem and the 
remarkable properties of the Levy stable probability laws serve as the basis of the LFs the- 
ory. Generalized central limit theorem says that the stable probability laws, like the Gaussian 
law, attract the distributions of sums of random variables [J. Due to this reason, Levy stable 
distributions naturally appear when evolution of a system, or result of an experiment are de- 
termined by the sum of independent, identically distributed random factors. The probability 
density functions (PDFs) of the stable laws exhibit slowly decaying, power-law asymptotic be- 
havior of the form , where a is called the Levy index, < a < 2. Due to this reason 
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the stable PDFs appear naturally in the description of random processes with large outliers, far 
from equilibrium. Another important property of LFs is their statistical self-similarity, or self- 
affinity 5 . Therefore, like the Brownian motion, they are naturally suited for the description 
of random fractal processes. LFs are ubiquitous in nature and they were observed in various 
fields of science, including physics (stochastic dynamics, turbulent flows [6171819) . turbulence 
and turbulent transport in magnetized plasmas |10lllll2fTB] . optical lattices [Tt], molecular 
collisions [15]), biology (heartbeat dynamics [16], firing of neural networks [17], searching on a 
folding polymer [IS], foraging movement [T^), seismology (recording of seismic activity [20]), 
stochastic climate dynamics [21], engineering (signal processing [22123124] ). economics (financial 
time series [25126127] V and even the spreading of diseases and dispersal of banknotes f28l. The 
list above is far from being complete; we only mention that recently illuminating experiments 
with LFs of light have been reported [29130) . 

The behavior of LFs in external fields brings a few surprising effects. The PDF in a harmonic 
potential evolves to a stationary state given by a Levy stable distribution 31,32. with the same 
stability index a as the underlying noise. For unimodal power law potentials steeper than a 
harmonic one U{x) oc |a;|'^, c > 2, the stationary state is characterized by the asymptotic power 
law decay, p{x) oc |a;|^^'^+"^^). These stationary probability densities decay faster [3334 35,36] 
than the corresponding Levy stable density. Moreover, the stationary PDFs in such potentials 
are bimodal, i.e. they possess a local minimum at the origin and two maxima at x ^ 0. The 
study of unimodal to multimodal bifurcations in this system during relaxation was addressed 
in Ref. [37 38 39 . For subharmonic potential c < 2 the decay p{x) slower than 

the decay of the corresponding Levy stable density. Moreover, stationary states only exist for 
c> 2-a 

The properties of the LFs in external fields were studied with the use of space-fractional 
Fokker-Planck equation and the Langevin equation for a particle driven by Levy noise. Mo- 
tivated by the model from stochastic climate dynamics [21[ . the barrier crossing problem for 
LFs in generic types of the potentials was studied in Refs. [41'42"43 I44I45I4()] . In particular, the 
method of decomposition of the Levy process into bounded jump component and compound 
Poisson part proposed in Refs. [45l46j allowed for treating the escape problem analytically (at 
least in the weak noise asymptotics), while solving the corresponding space- fractional Fokker- 
Planck equation poses considerable difficulties [12 • Further, the idea of decomposition was 
applied to describe a metastable behavior of a multi-well dynamical system driven by Levy 
noise [l^ . This have led to a new approach to fast simulated annealing problem [55] of efficient 
non-local search of the deepest potential well [49l50j . We will discuss the decomposition method 
in more detail below, see Sec. 13.11 

Levy processes lead to a richer behavior than Gaussian processes to which they reduce 
for a = 2. Moreover the limiting distributions for sums of independent identically distributed 
random variables do not need to be symmetric. The asymmetry of the noise can induce pre- 
ferred direction of the motion [51152) and asymmetry of stationary states [53) . This in turn is 
responsible for occurrence of the dynamical hysteresis |54l . 

In this paper we study, both analytically and numerically, LFs in a ratchet potential, i.e. 
in a periodic potential lacking reflection symmetry. Such spatial symmetry breaking gives rise 
to a net directed transport in presence of non-equilibrium fluctuations. The motion in ratchet 
potentials attracted considerable attention due to its role in fluctuation-driven transport, see 
for example Ref. [55] and references therein. By now, the physics of "Brownian motors" is well 
understood and the theory is well developed. However, not too much is known about Levy 
ratchets. In Ref. '56[ the authors solved numerically the Langevin equation with a white Levy 
noise for an overdamped particle in a ratchet potential. They demonstrated the appearance 
of directional transport using such measures of directionality as the position of the median 
of particle's displacements distribution characterizing the group velocity, and the interquantile 
distance giving the distributions' width. Such an approach allows the authors to study Levy 
ratchet in the whole range of the Levy index a. In Ref. [51] the analysis was restricted to 
the region 1 < a < 2, at the same time the numerical analysis of the Langevin equation is 
complemented by numerical solution of space-fractional Fokker-Planck equation. 
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The main goal of the present paper is to develop analytical theory of the Levy ratchet and 
to verify its predictions by numerical simulations. We have found that the probabilistic method 
of decomposition successfully applied to the escape problem for LFs in Refs. ,45^461 also works 
for the Levy ratchet description. As potential applications we can mention, following Ref. jSlj , 
a ratchet-like transport of impurities in magnetically confined fusion plasmas |57j , and a ratchet 
transport in biology, since biological systems are intrinsically far from equilibrium |58| . 

The rest of the paper is organized as follows. In Section [2] we formulate the Levy ratchet 
model. In Section[3]we give the essentials of the probabilistic decomposition method and present 
the results of analytical theory. In Section 2] we compare analytical results with the results of 
numerical simulations. The conclusions and summary are presented in Section [5] 



2 Model of the Levy ratchet 




Fig. 1. Exemplary potential used for inspection of the ratchet problem. The length of the potential seg- 
ment is controlled by L while the asymmetry of the potential is controlled by the potential asymmetry 
parameter g (0 < g < 1). For q = 1/2 the potential is symmetric. 



We study the following Langevin equation for an overdamped particle driven by the Levy 
stable noise L = dL/dt 

±{t)^-V'{x)+eL{t), (1) 

where V{x) is the potential, the prime denotes spatial derivative, and e is the amplitude of the 
noise, which will be considered to be small. The potential V{x) is periodic, V{x) = V{x + £), 
with potential wells flj, j = 0, ±1, ±2, .... The j-th potential well f2j is located at a; e (— -Cq -|- 
jC, C{1 — q) + jjC) where < q < 1, see Fig. [TJ As an exemplary potential V{x) we use 



V{x) 
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where L — a\ + a2, ai = (1 — g)£ and a-i — qC. The potential lacks reflection symmetry 
when q ^ 1/2. Compared with the Vq [sin(27ra:/£) -I- sin(47ra;/£)/4] potential typically used in 
literature, e.g. j55_, Eq. ^ offers an advantage of an easier control of the spatial asymmetry. In 
all the calculations presented below we use Vb = = 1. Figure [T] presents exemplary potential 
profile given by Eq. ([2]). 

In the integral form, Eq. ^ reads 



x{t)=x{0)~ I V'{x{s))ds + eL{t), (3) 

where L{t) is the Levy stable motion, that is the (formal) integral over the Levy stable noise 
in time. For a 1, the Fourier-transform of L{t), see Eqs. ^ and is 
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where a G (0,2] is the stability index, /3 e [—1,1] is the asymmetry (skewness) parameter, 
while c is a positive constant, c^/" being called a scale parameter. Initially, a particle starts 
its motion in the minimum of the potential, e.g. in the 0-th potential well. Furthermore, we 

assume (3 ^ ±1. 



3 Analytical approach to the Levy ratchet 

Analytical approach to the transport in a ratchet potential under action of Levy noise is based 
on the decomposition method, for which the main results were presented in |45| . whereas the 
mathematical details including proofs of corresponding theorems were given in |46j . Then, the 
dynamics in multi-well potentials was considered in |47| . Our exposition of the Levy ratchet is 
based essentially on these three papers. In Section [3] we first give an "intuitive" explanation 
of the decomposition method and then present theoretical results for main characteristics of 
ratchet transport, namely the mean exit time from a single well, a transition probability, a 
mean displacement and the current. 



3.1 Method of decomposition and single-well dynamics in case of symmetric Levy noise 

Let us consider the discretized version of Eq. ^ with the time step of integration At, 

x{n + 1) - x{n) = -V'{x{n))At + £,{n), (5) 

where ^(n) — L{nAt) — L{{n — 1)^^) — e{cAty/"(^{n) is the value of the noise variable, or the 
Levy jump, at the n-th interval of integration. Here, C{n) is the time-discrete Levy noise whose 
characteristic function is given by Eq. (U) with t — 1, c — 1. Thus, the PDF of these Levy 
jumps p = p{^) is characterized by the scale parameter, or characteristic length w — e{cAty^°' 
depending on the parameters of the noise and on the length At of the time interval. 

For the further exposition it is necessary to recall the asymptotics of the tails of the PDF 
p{£_) as ^ ±oo given e.g. in [SSJ Chapter 4.3]: 

PiO - --^^^ %eos(4^|r(^a)| gI^ - ^ ^ " 

Let us fix some level 6 > (depending on e) and decompose the noise ^{n) into two parts: 
Thus all ^(n) smaller in absolute value than 6 are considered to belong to the background part 
of the noise (containing most of C("-)), ^'iid large spikes with |^(n)| > S form a shot-like noise 
containing rare events, or outliers, see Fig. [2] showing the decomposition schematically. Since 
the values £,{n) are independent, the probability that a spike occurs on the step n does not 

depend on n and equals to Pgpiko = 

LooP(Od^ + I^PiOd^- For £ and At smah enough the 
characteristic length w is small, thus the probability Pspiko is determined by the large argument 
asymptotics of the Levy stable density given by Eq. ©, p{^) ~ A±w°'\^\-^^+°'\ ^ ±00, 
and behaves as Pspiko — {A/a)cAt£°'S~°', where A = A_ + A^. Therefore, the spikes form a 
Poissonian sequence of events, and, for 6 large enough, the spikes are well separated in time. 
The probability to have at least one spike on the unit time interval is 

P = —e'^S-". (7) 
a 

The probability not to have a spike on a time interval [0, t] is around {l — {Ac/a)e°'S~"Aty^^* ~ 
exp{-t/T) with 

r==-^e-"(5" (8) 
Ac 

being the mean time between two subsequent spikes. On the other hand, the background part 
of the noise has finite variance cr^ and its action over time intervals much larger than At can be 
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modelled by the action of a white Gaussian noise. The variance of the background component 
is given by ct^ = var(^(n)| ^ S) and, for large enough S, is again determined by the 

power-law asymptotics of the Levy stable density in its far tail: cr^ ~ const x e"(5^^"Z\t. The 
variance introduced by the background noise over the unit time interval (the background noise 
strength) due to the role of summation of variances for the white noise behaves as 

cr^ ~ const X e"(52-". (9) 

For e — !> it is always possible to chose 6 in such a way that both P in Eq. ((T]) and in 
Eq. ^ tend to zero. In Refs. [45146] the value of (5 = s^^^ was taken effectively to obtain the 
necessary estimates. 

Let us now fix small e and take S as discussed above. Between the two subsequent large 
spikes the particle is subjected only to a background noise, which, for small e, is so weak that 
the motion of a particle is almost deterministic, namely the particle slides towards the bottom 
of the potential well. The overall structure of such motion is well seen in the bottom panel 
of Fig. m For e small enough the background noise can be neglected, which would correspond 
to smoothing the rugged curve in the bottom panel of the figure and approximating it by a 
deterministic trajectory. On the other hand, the the time lag between the two subsequent spikes 
is so large that the particle reaches the bottom of the potential well between the two spikes. The 
only process of escape possible in this case, corresponds to the possibility that a strong enough 
spike "kicks" the particle from the bottom of the potential well, and an escape occurs in a single 
jump if its length exceeds the distance between the potential's minimum and maximum, see 
Fig. [5J The probability of jumping out of the well f^o is the conditional probability of having a 
positive spike larger than (1 — q)C or having a negative spike larger than qC in absolute value, 
that is 

P,.,, ^ Prob{e(n) < -qC or > (1 - q)C > s} ^ ^'^J ' ^ > 

where the wings of the Levy stable PDF representing the action of spikes per unit time are 
given by p(^) ~ A±w"^~^~°' , ±oo. Thus, we get 

/A_ A+ \ 6°' 

Pexit^ ^— + (11) 

which reproduces the result of Ref. [17] . 

Now we are able to calculate the mean exit time from a single well. Indeed, exit on a fc-th 
spike means that the first (k — 1) attempts were unsuccessful, ant the fc-th jump was big enough. 
The probability of this event approximately equals (1 — -Poxit)'^~^-Fcxit- The mean exit time in 
this case equals kT where the mean interspike time T is known, see Eq. ([8]). This yields 



(e)) ^ > fc"r(l - Pexit)''"'Pcxit - 



/ ^ '"-^ \^ ' exit; exit , y 

J T -'exit 
k—1 



2acos^|P(-a)| V(l - g)"" 9 



l + P 1-/3 



(12) 



X 



The picture of escape at small e differs drastically from the Kramers picture, namely instead 
of climbing up in the potential well the particle is thrown out from the well by a single large 
kick. Due to the fact that potential is translationally invariant the mean exit times (t*(£)) for 
all the potential wells are the same and equal to (T(e)). 

Now we can derive the asymptotic value of the transition probability Pi^i+k for a particle 
to make a jump from the i-th potential well to the (z -f fc)-th well. Assume for brevity that 
the particle starts at the well i7o, i = 0, and jumps to the well Qk, k ^ 1. The probability 
Pk — Po,k is then obtained similarly to Eq. ([T^ . Indeed, the transition to the well fit occurs if 
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the kick ^{n) overcomes the distance (1 — q)C + {k — 1)C but is smaller than (1 — q)C + kC. 
The probability of this event equals to 

Po^k ^ ProbjCH G [(fc - q)C, {k + 1- q)C] |e(n)| > s] ^ '^'"'^^ 



IlLpiOd^ + irPiOd^ (13) 
1 1 \A+ 3°" 



.{k-q)°' (k + l-q)"/ A ' 
Then with the help of the formula of the total probability we get 

n ~Vn p ^''-^p -P^^^ -n , ..^ (fc-g)-"-(fc + i-g)-" 

P.^^^ -P^-^ 1.(1 - P...) Po^. - - (1 + ^)(l-/3),- + (l + /3)(l-,)-- 

(14) 

Analogously, considering negative jumps we calculate transition probabilities pk — Pi,i+k for 
fc ^ — 1 as 

(1 - + + (15) 



3.2 Characteristics of transport 



A particle inserted into the i-th potential well spends a random time within the potential well, 
until it is kicked out of the well by a large spike, as it is described in subsection 13.11 It can be 
shown that in the small noise limit £ — > 0, the first exit time t^{s) from i-th potential fii is an 
exponentially distributed random variable with the mean value (T*(e)) = (T(e)), see Eq. (IT2]) . 
The transition probability Pi^i+k for a particle to make a jump from the i-th potential well to 
the {i + fc)-th potential well is given by Eqs. (fn]) and (|T5]) . 
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3.2.1 Splitting probabilities 

In order to elucidate the asymmetry of the first escape we use the splitting probability tt which 
is a probability of a first escape from a potential well to the right 

nn{a,P,q) = l^Pu= (i _ ^),-« + (i + _ ,)^. ' (16) 

A;^ 1 



or to the left 

7r^(a,/3,g) = )^ p. ^ (1 _ + (1 + /3)(1 - ' ^'^^ 



fc<-i 



Due to the asymmetry of the noise it is possible to find such a set of parameters for which the 
two values of splitting probabilities are equal. The balance condition 



^L(«,/3,g) =7rK(«,/3,g) = i (18) 



leads to the skewness parameter 



(1 - g)" + 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Q Q 

Fig. 3. The left panel presents values of the noise asymmetry leading to equal values of splitting 
probabilities, see Eq. (I19p . The right panel presents values of the noise asymmetry /3currcnt leading to 
zero current, see Eq. ((24}. 



3.2.2 Mean displacement 

For a > 1, the mean displacement of the particle at the time t is 

{X'it)) - {n)A, (20) 

where (n) — t/(T{e)) is the average number of jumps during time t. The mean value of the 
displacement A can be expressed as 

oc 

A = C J2 kpk^W+{a,l3,q)~W-{a,(5,q), (21) 

k— — oo,k^O 
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where 
and 

(l-/?)C(a,g) 



In above equations i^{a,q) denotes the Hurwitz zeta function, C(a, g) = X)feLo(^ ^" The 
mean jump length A, see Eq. (pOj) . and consequently the mean displacement {X'^{t)), see 
Eq. ipij) . are finite for a > 1 only. 

The main characteristics of the particles' transport in the Levy ratchet is the particle current. 
For the stability index a > 1, the current is defined as the time derivative of {X'^{t)), 

j{a,(3,q) = —{X it)} = — — . 

dt (T(e)) 

Consequently, for small noise intensities e — we get 

j(a,/3,(Z)cx£". (22) 



The multiplicative constant in Eq. ((22|) depends on parameters a and /? of the noise and on the 
asymmetry parameter q of the potential. 

We also determine the skewness parameter [3 of the noise, which leads to the zero particles' 
current. Indeed, the balance equation 

j(a,/3,q)=0 (23) 

yields the unique solution 

Pcurrcnt a, g) = 77 7— 77 T. 24) 



4 Numerical validation of the theory 

In order to show validity of the developed theory extensive numerical simulations have been 
performed for the stability index a € {0.7,0.9,1.1,1.5,1.9}, the noise asymmetry parameter 
13 G {0.0, ±0.5, ±0.9} and the potential asymmetry parameter q e {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8} 
with various (decreasing) noise intensities s. Numerical results were constructed by standard 
methods of integration of stochastic differential equations with respect to a-stable noises, see 
Refs. |60l61l62l63j . All numerical simulations have been performed with the time step of inte- 
gration At = 1Q~^, the number of realizations N = 10^, the potential depth Vq — 1, the scale 
parameter c = 1 and the potential segment length C = 1 . In general, all simulations in the limit 
of weak noise (small e) converge to theoretical predictions. From the whole set of simulations 
we have chosen exemplary results presented in Figs. H HllI 

Figure |4] presents the complementary cumulative distribution of the exit time 1 — CDF(t) 
(left panel) and the transition probabilities pk (right panel) for the stability index a = 0.9, 
the noise asymmetry (3 = —0.9 and the potential asymmetry q = 0.7. The left panel of Fig. U] 
demonstrates the exponential character of the exit time distribution whereas the right panel 
of Fig. U compares theoretical transition probabilities from the initial potential well SIq to the 
well i7fc given by Eqs. (fH)) and (jl5p (solid lines) with their numerical estimators (symbols), 
demonstrating perfect agreement between theoretical predictions and numerical simulations. 

Fig.[5]demonstrates the ratio of the numerically estimated values of the splitting probabilities 
7rf and their theoretical values tt]^ given by Eq. (fT7)) for different values of a. In the limit 
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500 1000 1500 2000 2500 1 10 100 



Fig. 4. The left panel demonstrates the exponential character of the (first) exit time distribution. The 
solid line in the left panel presents exp[— r/(r)^]. The right panel presents transition probabilities 
from the initial potential well (0) to the final potential well (fc). Circles and squares represent simulation 
results, while thin solid lines correspond to analytical formulas given by Eqs. (|14[1 and (|15|l . Simulation 
parameters: the stability index a — 0.9, the noise asymmetry /? — —0.9, the potential asymmetry 
q = 0.7 and the noise intensity e = 2~^. 




10'* 10'^ 10'^ 10"' 10° 

e 

Fig. 5. Ratio of simulated and theoretical values of splitting probabilities 7rf /tt^^. Simulation param- 
eters: the noise asymmetry /3 = —0.9, the potential asymmetry q — 0.7. 



of the small noise intensity the ratio between estimated and theoretical values of splitting 
probabilities tends to 1 indicating the agreement between the theory and simulations. For 
larger noise intensities numerical results deviate from predictions for weak noise regime. The 
sign of this deviation is different for a > 1 and for a < 1 . 

The left panel of Fig. |6] displays the ratio between the numerically estimated mean exit time 
{T{e))^ and the theoretical prediction (T(e)}^, see Eq. p^ . In the limit of small noise intensities 
e the ratio between estimated and theoretical values of mean exit times tends to 1 indicating 
that simulations performed corroborate theoretical findings. The right panel of Fig. [6] shows the 
scaling of the estimated mean exit time (T(e)} (symbols) with the noise intensity e along with 
the corresponding theoretical prediction, see Eq. p^ . shown as solid lines. Deviations from the 
predicted scaling visible for large noise intensities vanish for small values of e. 

Fig. [7] displays the splitting probability tt^ given by Eq. as a function of potential 
asymmetry q for the case of symmetric Levy noise with ^ = and different a. Here again the 
close agreement between analytical formulas and numerical results is observed. For symmetric 
noise and symmetric potential (q = 0.5) splitting probability is equal to 0.5. This indicates 
absence of the current in the system at hand. 
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Fig. 6. The left panel demonstrates the ratio of simulated and theoretical values of the mean escape 
time (t)^ / (t)'^ . The right panel presents the power-law scaling of the mean escape time as a function 
of the noise intensity e. Points represent simulation results while solid lines indicate the theoretical 
scaling of the mean exit time (r) with the noise intensity e. Simulation parameters: the noise asymmetry 
/3 — —0.9, the potential asymmetry q — 0.7. 



The left panel of Fig. [8] presents the mean displacement A defined for a > 1 only and given 
by Eq. (I21[) . and its numerical value estimated as 



(25) 



Ai being the number of the potential well to which a random walker was thrown from the initial 
one. For all values of a except for a = 1.1, which is too close to the boundary of convergence, 
the agreement between the numerical and the theoretical values is good. The right panel of 
Fig. |5] presents the splitting probability ttr as a function of the potential asymmetry q for the 
case of asymmetric Levy noise with /3 = 0.9. The numerical values (symbols) agree well with 
theoretical predictions (lines). 
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Fig. 7. Theoretical (solid lines) and simulated (points) values of the splitting probability -er as a 
function of the potential asymmetry q. Various curves correspond to different values of the stability 
index a. Simulation parameters: the noise asymmetry /? = 0.0, the noise intensity e — and the 
potential asymmetry q = 0.7. 



Fig. ini demonstrates the value of the splitting probability tt^ observed for the value of 
given by Eq. (jl9p (left panel) and values of the noise asymmetry parameter leading to 
equal values of splitting probabilities, see Eq. (IT9|) . (right panel). Results presented in Fig. [9] 
correspond to e = 2~^, which is larger than e in remaining figures. This is due to the fact that 
for noise asymmetry ji-^ given by Eq. (|19p longer simulation time is necessary. 
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Fig. 8. Theoretical (solid lines) and simulated (points) values of the mean displacement A'^ — (Ai) 
(left panel) and theoretical (solid lines) and simulated (points) values of the splitting probability ttr 
(right panel). Simulation parameters: the noise asymmetry /? — 0.9, the noise intensity s — and 
the potential asymmetry q = 0.7. 




Fig. 9. The left panel presents the simulated spitting probability ttl for the set of parameters for 
which the balance equality ttl = -kr — 0.5 should hold. The right panel presents values of the noise 
asymmetry leading to equal values of splitting probabilities, see Eq. ()19p . Simulation parameters: 
e = 2-«. 



The left panels of Figs. [10] and [TT] compare theoretical values of the current given by Eq. ([22]) 
with their numerical estimates as functions of the noise intensity e at time t — 2000. The power- 
law dependence j oc of the current, see Eq. (|22p is clearly demonstrated. For symmetric 
noises, see Fig. [TU] the absolute values of the current for q — 0.3 and q — 0.7 are the same. 
For q = 0.5, the theoretical value of the current vanishes, while numerical results fluctuate at 
a very low level. The right panels of Figs. [10] and [TT] present the dependence of the current j 
on time t for symmetric (/3 = 0) and asymmetric (/3 = 0.9) noise, respectively. Various curves 
correspond to different values of the potential asymmetry q. After a short transient the current 
reaches a constant value corresponding to a stationary regime of operation. 



5 Summary 

In this paper we have considered the transport properties of a particle embedded in a ratchet 
potential and subjected to a white Levy noise. In the case of weak noise the decomposition 
method proposed in Refs. |45l46j allowed us to calculate such characteristics of this trans- 
port as splitting probabilities of the first escape from a single well of a ratchet potential, the 
transition probability, the mean displacement of a particle and the particle current. In order 
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Fig. 10. The left panel presents the current j as a function of the noise intensity e at t = 2000 and 
various values of the potential asymmetry. The solid line presents theoretical scaling of the current on 
e. The right panel represents the current j as a function of time for the symmetric Levy noise and 
various values of potential asymmetry q. Simulation parameters: a — 1.5, 13 = 0, and e = 2~^. 
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Fig. 11. The same as in Fig. [10] for asymmetric Levy noise with /3 = 0.9. Simulation parameters: 
a = 1.5, e = 2"^ 



to confirm analytical results, we performed an extensive numerical simulations demonstrating 
good agreement with the theory at small values of the noise intensity. 

The authors acknowledge support by DFG within SFB555. AVC acknowledges financial support from 
European Commission via MC IIF, grant 219966 LeFrac. Computer simulations have been performed 
at Institute of Physics, Jagellonian University and Academic Computer Center, Cyfronet AGH. 
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